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We predict a resonance impurity state generated by the substitution of one Cu atom with a 
nonmagnetic atom, such as Zn, in the pseudogap state of a high-T c superconductor. The precise 
microscopic origin of the pseudogap is not important for this state to be formed, in particular this 
resonance will be present even in the absence of superconducting fluctuations in the normal state. 
In the presence of superconducting fluctuations, we predict the existence of a counterpart impurity 
peak on a symmetric bias. The nature of impurity resonance is similar to the previously studied 
resonance in the d-wave superconducting state. 



The effects of a single magnetic and nonmagnetic im- 
purity in high temperature superconductors have been 
studied intensively both theoretically and more re- 
cently experimentally by scanning tunneling microscopy 
(STM) Understanding of the impurity states in 

high-Tc materials is important because impurity atoms 
qualitatively modify the superconducting properties, and 
these impurity induced changes can be used to identify 
the nature of the pairing state in superconductors. 

Up to date theoretical analysis of the impurity states 
has been focused on the low temperature regime T <C T c 
well below the superconducting transition temperature 
T c . On the other hand it is well known that in the nor- 
mal state (T > T c ) of underdoped cuprates, the electronic 
states at the Fermi energy are depleted due to pseudo- 
gap (PG) Apg, as was seen by STM || and by angu- 
lar resolved photoemission ]To|j . One can consider the 
temperature evolution of the impurity state as the tem- 
perature increases and eventually becomes larger than 
T c . There are two possibilities for the evolution of im- 
purity resonance at T > T c : a) the impurity resonance 
gradually broadens until the superconducting gap van- 
ishes at which point the impurity resonance totally disap- 
pears and b) the resonance gets broader however survives 
above T c . Which of the possibilities is realized depends 
on the normal state phase into which the superconduc- 
tor evolves. It has been argued |]ll]jL2|] that in the un- 
derdoped regime the superconducting gap opens up in 
addition to the pseudogap present well above T c . We 
find that the impurity resonance survives above T c in the 
pseudogap state of high-T c materials. The position and 
the width of the resonance are determined by the impu- 
rity scattering strength and PG scale. In the absence of 
a PG above T c the impurity state disappears. 

The origin of PG state is one of the most strongly 
debated issues. Some models attribute PG to supercon- 
ducting phase fluctuations above T c |l3j ; others to a com- 
peting non-superconducting order parameter p[ . An- 
other possibility is that even within the PG regime there 



are at least two distinct sub-regimes - strong and weak 
pseudogaps, with weak pseudogap occurring at higher 
temperatures due to antiferromagnetic fluctuations, and 
strong pseudogap being related to superconducting fluc- 
tuations |^@. 

In this article we address the impurity induced reso- 
nance or quasibound state that is generated by a strong 
nonmagnetic impurity scattering in a CuO plane in the 
normal state of high-T c materials. Specifically we calcu- 
late the resonant state generated by the substitution of 
one Cu atom with a Zn atom using the self consistent 
T matrix approach. We rely on the fact that the den- 
sity of states (DOS) is depleted at the Fermi energy in 
the PG regime. We argue that the mere fact that DOS 
is depleted at the Fermi energy is sufficient to produce 
impurity resonance near the nonmagnetic impurity, such 
as Zn. However no particular use of the superconducting 
correlations above T c has been made in our analysis. For 
example, the results we present, will be valid in the PG 
state with no superconducting phase or amplitude fluc- 
tuations above T c , as long as there are interactions that 
lead to the PG state, as indicated by a depleted DOS. 
This is an important caveat that broadens the validity of 
the model regardless of the microscopic origin of the PG 
in the high-T c superconductor. The approach we take is 
similar to the previous analysis of the nonmagnetic im- 
purity in the superconducting state JlJ]. See also figure 

The superconducting fluctuations are not required for 
the formation of the impurity state in the PG regime. 
However, in the presence of superconducting fluctuations 
an additional important feature of the impurity state is 
expected to appear. The natural quasiparticles in the 
superconducting state are the Bogoliubov quasiparticles, 
which are a linear combination of a particle and a hole. 
As a consequence, an impurity state in a superconductor 
appears both on positive (particle) and negative (hole) 
biases p]Jl7j • In the phase-fluctuating pseudogap the Bo- 
goliubov quasiparticles acquire a finite lifetime. However, 
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the particle-hole symmetry of the impurity states should 
remain to the extent to which the superconducting quasi- 
particles are defined [0. Hence, the PG impurity states 
can serve as an extremely local probe able to distinguish 
between the superconducting and non-superconducting 
scenarios for the PG regime. 
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FIG. 1. An impurity state in a high T c superconductor: 
(a) The DOS in the pseudogap regime used in this article 
(see also [11]) and (b) the DOS in the superconducting state 
as was used in [1]. In both phases there is a resonant state. 



To be specific we need a model DOS that captures the 
main features of the PG in high-T c materials. For this 
purpose we use the DOS that was measured Q by Lo- 
ram et al.. In this work it has been argued that the DOS 
is a linearly vanishing function of energy within the Apq 
energy range near the Fermi surface, see figure 1(a). 

We find that such a model indeed gives rise to an impu- 
rity bound state with energy f2' and decay rate Cl" equal 
to 
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where we have assumed the impurity scattering to be 
strong enough so that the result can be calculated to 
logarithmic accuracy with In |2C/A^o | > 1 [Q. This is 
the main result of our work, which we will derive in the 
remainder of this paper. 

The Hamiltonian for this single impurity problem is 
given by 



Uh 



kk' cr 



Ck'a, 



(2) 



where U is the strength of the scalar impurity potential 
taken to be nonzero only at r = 0. The scattering matrix 
M can be written as 



T = 



U 



u 



(3) 



G (Q) 



(4) 



with Go(w) the on site Green's function [[19]]. The states 
generated by the impurity are given by the poles of the 
T matrix: 

1 

U' 

This is an implicit equation for Cl as a function of U, the 
strength of the scattering. This solution can be complex, 
indicating the resonant nature of the virtual state. To 
solve this equation, we split Go into its imaginary and 
real part G = G + iG '. But also G '(w) = -ttN (lj) 
with Nq(ui) the density of states. 

Measurements on the electronic specific heat by Loram 
et al. (jllj show that the normal state pseudo gap opens 
abruptly in the underdoped region below a hole doping 
equal to p cr i t ~ 0.19 holes/Gu02- Inspired by these data, 
we will assume that around the pseudogap region, states 
are partly depleted and the density of states is linear, that 
is N(u) = N \lu\/A pg for \u\ < A PG and N(lu) = N for 
Ap G < \w\ < W/2 with W the bandwidth. This density 
of states is depicted in figure ||(a) . As it is obvious from 
the solution of equation (|4|) , the precise position and the 
width of the resonance will depend on the specific form 
of the PG. We will use this linearly vanishing PG DOS. 
Results for the other form of N (oj) can be obtained in a 
similar fashion |po| . 

From the Kramer-Kronig relation |2^] 
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with P Cauchy's principle value, one can calculate the 
real part G' giving 
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This function is depicted in figure y(b) together with 
l/U. If 2UNo > e, with e Euler's constant, one can see 
from this figure that equation (jij) has four solutions. But 
because the width of a resonance state is proportional 
to | Q |, the only state with sharp width is the solution 
with \Cl\ close to zero and we will only consider this so- 
lution. After expansion in ui of equation (JsJ) we arrive at 
an expression for this solution CI of equation (0): 
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which is solved to logarithmic accuracy by expression 
(0). Using this formula, taking A^ = 1 state/ eV, 
Ape ~ 300AT ~ 30meT^ and the scattering potential 
U ~ ieV, we estimate Ct ~ 20^ ~ 2meV as was found 
by Loram et al. (ll]] This energy is close to the Zn res- 
onance energy luq = 16K, seen in the superconducting 
state 0. 
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FIG. 2. (a) Fhe density of states N(w) = -G'o(u)/n. 
Around the pseudogap states are only partly depleted e.g. 
N(u) = No\uj\/Apg, where N(u) = No for A pg < M < W/2 
with W the bandwidth, (b) The real part G' (u) of Green's 
function together with 1/U and U positive. SI' is the real 
part of the solution of the equation Go (SI) = 1/U close to 
zero and therefore with sharp bandwidth, (c) The impurity 
induced resonance at Si' — — Apg/2Z7JVo ln(2t7jVo). Because 
the other three solutions of equation (Q) have much broader 
bandwidth, they are not depicted here. All the figures are 
taken on the impurity site. 



. This calculation would require a specific model for 
the PG state and goes beyond the scope of this paper. 

While no superconductivity is required to form the im- 
purity state in the PG, if the superconducting fluctua- 
tions are present then an additional satellite peak should 
appear on a symmetric bias due to the particle-hole na- 
ture of the Bogoliubov quasiparticles. The relative mag- 
nitude of the particle and the hole parts of the impurity 
spectrum can be used to determine the extent to which 
the PG is governed by the superconducting fluctuations. 
In the case of fully non-superconducting PG there should 
be no observable counterpart state. An optimal impurity 
for such determination would appear to be Ni, which 
breaks the particle-hole symmetry more weakly than Zn 
even in the superconducting state. Combined with other 
experimental proposals 22 2j| , the impurity state can 



help to better understand the mysterious PG state. 

In conclusion, we find the resonance state that is in- 
duced by the nonmagnetic impurity in the normal state 
of a high-T c superconductor in the PG regime where the 
DOS at the Fermi surface is depleted. For the particular 
model of linear DOS we find that impurity state energy 
is given by equation (|l|) . Impurity states survive at high 
temperature T > T c since the PG produces the DOS de- 
pletion. This depletion is all that is necessary to produce 
the intragap state. In our solution we did not rely on a 
superconducting phase fluctuations above T c to generate 
the impurity state. We estimated the energy of the Zn 
impurity resonance to be at 20K, assuming a impurity 
potential equal to U = AeV . This resonance energy is 
indeed close to the energy of the Zn resonance at 16K 
in superconducting state 
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The solution of the impurity state deep in the super- 
conducting regime involves two aspects: the energy po- 
sition and the width of the resonance and secondly, the 
real space shape of the impurity state. We have discussed 
the energy of the impurity state above. Great advantage 
of the on-site impurity solution for the localized potential 
U is that only on-site propagator Gq(uj) enters into cal- 
culation. Hence the knowledge of the DOS was sufficient 
to calculate the impurity state. On the other hand, to 
calculate the real space image of impurity induced res- 
onance, one would require more detailed knowledge of 
the Green's functions in the PG regime. Quite gener- 
ally, one would expect for a d-wave like PG with nearly 
nodal points along the (±tt/2, ±tt/2) directions, that the 
impurity resonance in the pseudogap regime would be 
four-fold symmetric, similar to superconducting solutions 
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